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ABSTRACT 

Non-linear effects in the evolution of the axion field in the early Universe may lead 
to the formation of gravitationally bound clumps of axions, known as "miniclusters." 



! Minicluster masses and radii should be in the range Mine ~ 10-^^0 and R 

mc ~ 10 cm 



and in plausible early-Universe scenarios a significant fraction of the mass density of the 
^ ■ Universe may be in the form of axion miniclusters. If such axion miniclusters exist, they 

^ would have the physical properties requhed to be detected by "femtolensing." 
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Ph' Subject headings: dark matter — axions; gravitational lensing — femtolensing 
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■ 1. Introduction 



Gravitational lensing is now a well established field of astronomy, with a wide variety of 



' astrophysical and cosmological applications. Among the many interesting and important implications 
^ of lensing phenomena is the possibility of interference effects between lensed images of point-like 
astrophysical sources ([Mandzhos 198 1| ; Phanian 1983| ; Schneider fc Schmidt-Burgk 1985| ; peguchi"fc 



[Watson 19861 ; [Peterson fc Falk 1991|) . If the angular separation of the images of a point-like source 
is comparable to the wavelength of the light, interference between the lensed images will lead to a 
distinctive fringe pattern both in coordinate space and in energy spectrum. The interference pattern 
in the energy spectrum will depend upon the location of the detector. Since well separated detectors 
will see different patterns, there is a distinctive signature of the effect. 

If gamma-ray bursters are of cosmological origin they would be small enough to act as point-like 
sources, and lenses of extremely small masses, 10^^^ M© < M < 1O~^^M0,[] would produce a 
frequency-dependent signal in their spectrum ( [Gould 1992| ; [Stanek et al 1993[) . For source and lens at 



cosmological distances, the angular separation of the images would be in the femto-arcsec range, so the 
phenomenon was dubbed "femtolensing" by Gould (1992). 



"'^Post-WKB effects may extend the possible detection range of the lens up to 10 Mq (Ulmer & Goodman 1995). 
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For source and lens of cosmo logical distances, the Einstein ring radius is Re ~ 
3 X lO^°(M/lO~^^M0)^/^cm. This means that if two gamma-ray burst detectors are separated 
by distances larger than Re, typically one detector will be within the Einstein ring and the other will 
not. Measurement of different fluxes in two widely separated detectors would extend the detectable 
range of lens masses to M < 1O~^M0 ( [Nemirofl:" fc Gould 19951 ), which can be called picolensing. Thus, 
femtolensing and picolensing could be a probe of dark matter in a mass range difficult to probe by 
other means, not to mention settle the issue of the origin of gamma-ray bursters. 

Although the phenomenon of femtolensing (henceforth we omit picolensing for brevity, but 
everywhere below we mean both phenomena) is quite intriguing, the possibility has not attracted a 
great deal of attention because of a lack of attractive femtolens candidates. Among the front runners in 
an admittedly weak field are snowballs, black holes, and small molecular clouds. Snowballs (lumps of 
primordial H and He) are generally believed to evaporate on time scales shorter than the Hubble time 
( Pe Rujula et al 1992|) and for this reason they are not very attractive. A significant fraction of closure 
density in black holes in the requisite mass range seems difficult to arrange without a very peculiar 
spectrum of primordial fluctuations to prevent too many small holes which would be evaporating now. 
As an unattractive possibility in a field of even more unpalatable choices, small molecular clouds are 
usually considered the best bet for femtolensing. 

Even if one thought of a way to arrange baryons into femto-mass objects, it is unlikely that 
they would contribute a large enough fraction of closure density for a reasonable optical depth for 
femtolensing. Even if all baryonic matter would form femtolenses, their density is constrained by 
big-bang nucleosynthesis to be less than about 10% of closure density. 

Since non-baryonic dark matter is expected to dominate the mass density of the Universe, 
perhaps it is easier to arrange a near critical density of non-baryonic material into objects capable of 
femtolensing. But most particle dark matter candidates are very weakly interacting, and they would 
not be expected to clump on the femtolensing mass and density scales. 

In this paper we shall show that the situation may be different if axions are the dark matter. Here 
we outline a scenario where a significant cosmological density of axion miniclusters (|Hogan fc Kess 



1988n iKolb fc Tkachev 1993| , 1994a) form with the necessary properties to be natural candidates both 
for femto- and picolensing. 



2. Formation of Axion Miniclusters 



The invisible axion is among the best motivated candidates for cosmic dark matter. The axion 
( [Weinberg 1978| ; [Wilczek 1978|) is the pseudo-Nambu-Goldstone boson resulting from the spontaneous 
breaking of a U{1) global symmetry known as the Peccei-Quinn, or PQ, symmetry introduced to 
explain the apparent smallness of strong CP- violation in QCD ( |Peccei fc Quinn 1977|) . 



There are strin gent astrophysical (picus et al 1978|; [b'ukugita et al 1982 ; pdatsuda fc Yoshimura 



T988| ; [Turner 1988] [Raffelt fc Seckel 1988| ; |Mayle et al 1988|) and cosmological ( Preskill et al 1983 



[Abbott fc Sikivie 1983| ; [Dine fc Fischler 1983[) constraints on the properties of the axion. In 



particular, the combination of cosmological and astrophysical considerations restrict the axion decay 
constant fa and the axion mass to be in the narrow windows 10^° GeV < fa 10^^ GeV, and 
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10~^ eV < rria < 10^^ eV.0 The contribution to the mean density of the Universe from axions in this 
window is guaranteed to be cosmologically significant. Thus, if axions exist, they will be dynamically 
important in the present evolution of the Universe. 

The axion is effectively massless well above the QCD confinement temperature of a few hundred 
MeV. At the confinement scale and below, QCD effects produce a potential for the axion of the form 
V{d) = m^/f [1 — cos(6')], where the axion field was parametrised as a dimensionless angular variable 
6. Above the QCD scale there is no reason for the axion field to be at its low-temperature potential 
minimum, so as the axion potential turns on the field must relax to the minimum of the potential. Thus 
the energy density in axions corresponds to coherent scalar field oscillations, driven by a displacement 
of the initial value of the field (the "misalignment" angle) away from the eventual minimum of the 
temperature-dependent potential. 

There will be nothing special in the subsequent evolution of the axion field if the initial 
misalignment angle is constant throughout the Universe. This may in fact be a reasonable picture if the 
Universe underwent inflation with the "re-heat" temperature well below fa- However, equally plausible 
early-Universe scenarios lead to fluctuations of order unity in the misalignment angle on scales larger 
than the Hubble radius as the potential starts to develop. Since the final energy density depends upon 
the initial misalignment angle, spatial fluctuations in the misalignment angle are transformed into 
spatial fluctuations in the axion density, which later lead to tiny gravitationally bound "miniclusters." 

It is easy to understand that miniclusters will be very dense objects. Background axion density 
scales with temperature as Pa{T) = 3Teqs/4 for T <C Aqcd, where s oc is the entropy density and 
Teq is the temperature of equal matter and radiation energy densities. Consider a region with axion 
over-density, pa = {1 + ^)pa, where $ is larger than unity. The matter density in that region will 
dominate the radiation density at a temperature T<i, = (1 + $)Teq. At this time the density fluctuation 
became non-linear, separate out from the cosmological expansion, gravitationally collapse, relax, and 
form a minicluster with the approximate density it had at T$. A detailed study of this leads to a final 
minicluster density of ([Kolb fc Tkachev 1994b|) 

p^c ^ 140$3(1 + $)Pa(req) ^ 3 X 10-'^^''{1 + $) [nah^f g Cm'l (1) 

Even a relatively small increase in $ is important because the final density depends upon $^ for $ > 1. 

In Hogan & Rees (1988) it was assumed that typical values of $ corresponding to a minicluster are 
of order unity. The typical mass of a minicluster in Hogan & Rees (1988) was assumed to correspond 
to the mass of all axions inside the horizon at T ~ 100 MeV, with the result M^c ~ 10"^ M©. 



In our numerical investigations of the dynamics of the axion field around the QCD epoch ( |Kolb 



& Tk achev 1993| , 1994a) we found that as the oscillations commence, important, previously neglected. 



non-linear effects of the field self-interaction can result in the formation of transient soliton-like objects 
we called axitons. The non-linear effects result in regions with $ much larger than unity, possibly as 
large as several hundred, leading to enormous minicluster densities. We also found that the minicluster 
mass scale is set by the total mass in axions within the Hubble radius at a temperature around 
T = Ti 1 GeV when axion mass is equal to H. This lowers the minicluster mass from the estimate 
of Hogan and Rees to about 1O~^^M0. 



^ For a review see Turner 1990 or FLaffelt 199C 
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Our previous calculations suggested that miniclusters could be a candidate for femtolensing. We 
have since extended and improved our original calculation of the properties of axion miniclusters 
using the same approach as in [Kolb fc Tkachev 1993| , but with better numerical resolution and 
using a different spectrum of initial fluctuations of 9. The new calculation confirmed the basic 
picture of minicluster formation and sharpened the predictions of the masses, radii, and abundance of 
miniclusters. In the next section we summarize the results relevant for femtolensing. 



3. Physical Properties of Axion Miniclusters 

The evolution of the axion field in an expanding, spatially fiat Universe during the epoch when 
axion mass switches on is governed by a wave-like equation with a non-linear, time-dependent potential 
( IKolb fc Tkachev 1993| ): 

i)" - A^ij + r]"+^ sm{^/r]) = , (2) 

where ■?/' = V^y 4'" = (Pip/drf, t] is the conformal time (normalized to t] = 1 when the inverse of the 
axion mass is equal to the Hubble radius — i.e., i]{Ti) = 1), and A is the Laplacian with respect to 
comoving coordinates. The factor of n in this equation arises from the way the potential develops due 
to instanton effects: 'm'^{T]) = m'^{r] = 1)1]", with n = 7.4 ± 0.2 ([Gross et al 1981]) . 



As an initial conditions in the new calculation we considered the possibility that the main source 
of axions is provided by the decay of axion strings prior to T = Ti ( |Davis 198(j| ; [Harari fc Sikivi 1987| ; 



Battye fc Shellard 1994| ). This might be expected if the reheating temperature after inflation is larger 



than the temperature of the Peccei-Quinn phase transition. In such a scenario the spectral energy 
density of axions produced by string decay corresponds to a power law: dpa/du oc uj~^, where u = \k\ 
(recall that the axion effectively is massless prior to ?7 = 1). This corresponds to \a{k)\ = const x 
for the Fourier amplitudes of the axion field. In our latest calculations we used initial conditions at 
rj = 1 with amplitudes of the Fourier modes corresponding to the above spectrum and with random 
phases. Overall amplitude was chosen to ensure rms value of 9 to be ir/^/S. 

The final density distribution of axions in one of the 2-dimensional slices through the integration 
volume is shown in Figure 1. Three high amplide density peaks which are seen on this figure 
results from non-linear interaction and are absent in the case of the harmonic potential. They form 
miniclusters of interest at T ~ $Tcq. The evolution of the axion field with the above initial conditions 
was qualitatively similar to the evolution using the white-noise initial conditions of our previous 
calculations. From the results of the numerical calculation we determined the mass, radius, and 
abundance of axion miniclusters. 

3.1. Minicluster mass. The mass of axion minicluster is set by the scale when axion 
oscillations commence. As mentioned above, we have defined Ti as a temperature when the 
axion mass equals the Hubble constant, ma(Ti) = H{Ti). Using results of Turner (1986) we use 
Ti = (/„/10^2QgY) -0.175 (^^^^/200MeV)°-^ GeV. In what follows we assume Ti = 1 GeV. As a 
reference mass scale we define the mass of all axions in a cubic volume of length equal to the Hubble 
length at Ti, and density equal to the mean cosmological density: Mi = 'K'^g^.{Ti)TfTec^/?)QHl, where 
g^{Ti) ^ 62 is the number of an effectively massless degrees of freedom at temperature Ti. Using 
H = 1.66^y^rVMpi we find Mi = 8.2 x lO^^^^a/i^ M©. 

Our numerical integration of equation (0) shows that there is a spectrum of minicluster masses, 
but typically the mass of the axion minicluster is a reasonable fraction of Mi, Mmc ^ O.IMq. We 
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estimate the average minicluster mass to be of order M^c ~ 10"'^ Mq. 

Masses of miniclusters are relatively insensitive to the particular value of $ associated with the 
minicluster. 



3.2. Minicluster radius. Typical axion miniclusters have a large density contrast, $ > 1 ( |Kolb fc 



Tkacl ^ev" 1993|) , see Figure 1. Equation (|ID gives the minicluster radius as a function of M and $: 

1 /3 

^ 2 X 101^$-^ (1 + <|.)-^/^ {n^h'Y' {yo^m^) cm. (3) 

The radius of a particular minicluster depends upon its associated value of $. In the limit $ > 1 
the radius of a minicluster scales as 

3.3. Minicluster abundance. From the numerical calculation we have extracted the mass fraction 
of axions in density peaks of $ larger than a given value $o- It is plotted for several values of conformal 
time f] = 1,3, and 4 in Figure 2. By the time ?7 > 4 we are no longer able reliably to resolve the highest 
density peaks because density contrasts became too large for our 256^ grid. Although at ?7 = 4 the 
evolution has not completely frozen out, (especially at the large-<l> end) we can set a lower bound for 
the density fractions, e.g., more than 13% are in miniclusters with $ > 10 and 70% of all axionic dark 
matter are in miniclusters ($ > 1), see Figure 2. 



4. Miniclusters as Femtolens 

For axion miniclusters (or, for that matter, any other possible candidate) to be detected by 
femtolensing (or picolensing) the following three conditions must be satisfied: 



i.) The minicluster has to be in the mass range, 10 Mq < ^ 10 Mq ( |Gould 1992| ; |U Imer 



Goodman 1995 ) for femtolensing, or M^c ^ 10 ^Mq ( [Nemiroff fc Gould 1995|) for picolensing. 



a.) The physical radius of the minicluster has to be smaller than its Einstein ring radius. For 
source and lens of cosmological distance the Einstein ring radius is, see e.g. Press & Gunn (1973): 

where as usual Hq = lOOh km s~^Mpc~^. 

Hi.) The third condition is that a significant fraction of the total mass density of the Universe 
must be in the form of miniclusters, say fimc ~ 0.1, for the lensing rate to be reasonable. This condition 
is independent upon the lens mass ( Press fc Gunn 1973|) . 



These conditions are non-trivial, which is why there were no serious candidates for femtolensing. 
But we shall now demonstrate that all three conditions are naturally met for the case of axion 
miniclusters. 

That miniclusters satisfy the mass condition follows from the fact that the minicluster mass is 
determined by the total mass of axions within the Hubble radius at T ~ 1 GeV. Of course there will 
be some spread in minicluster masses, but as discussed in the previous section, the masses should be 
peaked around lO^^^M©, well within the desired range for lensing. 
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Using the radius of miniclusters from equation (|^) and the Einstein ring radius from equation 

(ID, we find RE/Rmc ^ 0.15$ (1 + <l>y^^ Qah^^^ {M/10-^^MQ f^\ The requirement that miniclusters 
be smaller than their Einstein ring radius places restrictions on $ as a function of M. We see that if 
fla ~ 1; then miniclusters with M ~ 10~^^Mq and $ > 5 — 10 (depending upon h) are within their 
Einstein ring radius. Actually, objects with the radius comparable to their Einstein radius have a 
better chance to be detected as femtolenses (compared to point mass objects) if they are in the mass 
range 10~^^Mq < M < 10~^^ Mq ( |Ulmer fc Goodman 1995| ). From Figure ^ we see that more than 
about 20% of all dark matter axions are in miniclusters with $ > 5 and will be within their ring radius. 

We conclude that if axions are the dark matter with Qa close to unity, then the fraction of Q in 
miniclusters capable of femtolensing is 15% - 20% (see Fig. H), i. e. will be in excess of 10% — large 
enough to result in a detectable rate. 

Since miniclusters in the IGM as well as in galactic halos contribute to femtolensing events, 
questions such as: how are they clustered, do miniclusters survive the epoch of galaxy formation, are 
they disrupted in their encounters with stars, etc., may be irrelevant here. Since large-'l' miniclusters 
are very dense, form early, and are well separated from each other, they should escape mutual 
disruption and merging. 

It can be important for the present discussion that miniclusters with $ > 30 relax to Bose-stars 
( P?kachev 1991| ; [Kolb fc Tkachev 1993| ) and consequently became even densier and compact. 



We conclude that axion miniclusters are candidates for femtolensing and picolensing. We do not 
believe that the conclusion is sensitive to specific initial conditions, but is generic so long as there are 
large super-horizon-sized fluctuations in the misalignment angle during the development of the axion 
potential. 



We thank S. Colombi, A. Gould and A. Stebbins for useful discussions. This work was supported 
by DOE and NASA grant NAG 5-2788 at Fermilab and DOE grant DE-AC02-76ER01545 at Ohio. 
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Fig. 1. — This is a two-dimensional slice through a three dimensional box showing the distribution of 
axion energy densities for T ^ Tgq. The height of the plot corresponds to $ = Spa/pa = 20, and the 
width to a length of 4if~^(Ti) (this corresponds to a comoving length of about 0.2 pc). 
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Fig. 2. — The mass fraction of axions in miniclusters with $ greater than $o as a function of $o- By 
T) — A the evolution has very nearly frozen. 
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